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Exploiting Known Interference as Green Signal 
Power for Downlink Beamforming Optimization 
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Abstract —We propose a data-aided transmit beamforming 
scheme for the multi-user multiple-input-single-output (MISO) 
downlink channel. While conventional beamforming schemes 
aim at the minimization of the transmit power subject to 
suppressing interference to guarantee quality of service (QoS) 
constraints, here we use the knowledge of both data and channel 
state information (CSI) at the transmitter to exploit, rather 
than suppress, constructive interference. More specifically, we 
design a new precoding scheme for the MISO downlink that 
minimizes the transmit power for generic phase shift keying 
(PSK) modulated signals. The proposed precoder reduces the 
transmit power compared to conventional schemes, by adapting 
the QoS constraints to accommodate constructive interference 
as a source of useful signal power. By exploiting the power of 
constructively interfering symbols, the proposed scheme achieves 
the required QoS at lower transmit power. We extend this 
concept to the signal to interference plus noise ratio (SINK) 
balancing problem, where higher SINK values compared to the 
conventional SINK balancing optimization are achieved for given 
transmit power budgets. In addition, we derive equivalent virtual 
multicast formulations for both optimizations, both of which 
provide insights of the optimal solution and facilitate the design of 
a more efficient solver. Finally, we propose a robust beamforming 
technique to deal with imperfect CSI, that also reduces the 
transmit power over conventional techniques, while guaranteeing 
the required QoS. Our simulation and analysis show significant 
power savings for small scale MISO downlink channels with 
the proposed data-aided optimization compared to conventional 
beamforming optimization. 

Index Terms —Downlink beamforming, convex optimization, 
power minimization, SINK balancing, constructive interference 

I. Introduction 

The power efficiency of wireless transmission links has 
recently attracted considerable research interest, in-line with 
the global initiative to reduce the CO 2 emissions and op¬ 
erational expenses of communication systems. Accordingly, 
transmit beamforming schemes for power minimization have 
become particularly relevant for the downlink channel. Capac¬ 
ity achieving non-linear dirty paper coding (DPC) techniques 

Manuscript received May 30, 2014; revised September 19, 2014 and 
February 9, 2015; accepted April 18, 2015. The associate editor coordinating 
the review of this manuscript and approving it for publication was Prof. Y. 
W. Hong. 

Christos Masouros is with Department of Electronic & Electrical Engi¬ 
neering - University College London, Torrington Place, London WC1E 7JE, 
UK, E-mail: chris.masouros@ieee.org. 

Gan Zheng is with School of Computer Science and Electronic 
Engineering, University of Essex, Colchester, C04 3SQ, UK, E-mail: 
ganzheng(a)essex.ac.uk. He is also affiliated with Interdisciplinary Centre 
for Security, Reliability and Trust (SnT), University of Luxembourg, Luxem¬ 
bourg. 

This work was supported by the Royal Academy of Engineering, UK and 
the Engineering and Physical Sciences Research Council (EPSRC) project 
EP/M014150/1. 


0 0 have been proposed for pre-subtracting interference 
prior to transmission. The DPC methods developed so far 
are in general complex as they require sophisticated sphere- 
search algorithms (3) to be employed at the transmitter and 
assume codewords with infinite length for the encoding of the 
data. Their suboptimal counterparts 0 © offer a complexity 
reduction at a comparable performance. Still however, the 
associated complexity is prohibitive for their deployment in 
current communication standards. On the other hand, linear 
precoding schemes based on channel inversion (7)- (TTJ of¬ 
fer the least complexity, but their performance is far from 
achieving the optimum maximum likelihood bound. Their non¬ 
linear adaptation, namely vector perturbation (VP) precoding 
G3-GD provides a performance improvement at the expense 
of an increased complexity since the search for the optimal 
perturbation vectors is an NP-hard problem, typically solved 
by complex sphere search algorithms at the transmitter. 

More relevant to this work, optimization based techniques 
that directly minimize the transmit power subject to quality of 
service (QoS) constraints - most commonly to interference- 
plus-noise ratio (SINR) - have been studied for broadcast 
channels in 03. where convex optimization problems of such 
nature were proposed. Recent works have focused on the 
utility maximization in MIMO interfering broadcast channels 
03.09 and full duplex radio (20). For the cases of channel 
state information (CSI) errors, robust versions of these opti¬ 
mizations have been studied in [21] - [30). In [25], a robust 
max-min approach was developed for a single-user MIMO 
system based on convex optimization. Later in fTT) , the robust 
transmission schemes to maximize the compound capacity for 
single and multiuser rank-one Ricean MIMO channels were 
addressed, based on the uncertainty set in [[26). Robust beam¬ 
forming for multiuser multiple-input single-output (MISO) 
downlink channels with individual QoS constraints under an 
imperfect channel covariance matrix was studied in (16), (27) . 
Recently in (28) , the optimal power allocation over fixed 
beamforming vectors was obtained in the presence of errors 
in CSI matrices. Most recently, efficient numerical solutions 
to find conservative robust beamforming for multiuser MISO 
systems with mean-square-error (MSE) and SINR constraints 
and different bounded CSI errors have been developed in 
(29) , (30) . Moreover, SINR balancing optimizations have been 
proposed in (D where the minimum achievable SINR is 
maximized, subject to a total transmit power constraint. 

This paper is based on the concept of interference exploita¬ 
tion, first introduced in © 0 (32), [33] where analytical in¬ 
terference classification criteria and low-complexity precoders 
based on channel inversion were derived. The analysis showed 
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how the knowledge of both CSI and data at the base station 
can be used to predict and classify interference into construc¬ 
tive and destructive. The closed-form precoders in ® © 
showed that by exploiting the constructive part of interference, 
higher receive SNRs can be provided without increasing the 
per user transmit power. However, when precoders are fully 
optimized, less is understood about the performance gain of 
the constructive interference approach over the conventional 
optimization, e.g., (TS). Accordingly, in this work we aim 
to improve such optimization-based precoders by exploiting 
constructive interference as a source of signal power. By 
doing so, we further reduce the transmit power required for 
guaranteeing the SNR constraints in the optimization, thus 
improving the power efficiency of transmission. For clarity, 
we list the contributions of this paper as follows: 

1) We introduce a new linear precoder design for PSK that 
a) reduces the transmit power for given performance 
compared to existing precoders based on the proposed 
constructive interference regions, and b) as opposed to 
conventional precoders, applies to scenarios with higher 
numbers of users than transmit antennas, 

2) We further adapt this concept to the SINR balancing 
problem, where higher SINR values compared to the 
conventional SINR balancing optimization are achieved 
for given transmit power budgets, 

3) We re-cast the optimization problem as a virtual multi¬ 
cast optimization problem based on which we derive the 
structure of the optimal solution and develop an efficient 
gradient projection algorithm to solve it, 

4) We use the multicast formulation to derive a robust pre¬ 
coding scheme suitable for imperfect CSI with bounded 
CSI errors. 

We note that, while the following analysis focuses on PSK 
modulation, the above concept and relevant optimizations can 
be extended to other modulation formats such as quadrature 
amplitude modulation (QAM) by adapting the decision thresh¬ 
olds of the constellation to accommodate for constructive 
interference (32j. It should be stressed however, that the pro¬ 
posed schemes are most useful in high interference scenarios 
where low order modulation such as BPSK and QPSK are 
used in the communication standards to ensure reliability (34| . 
In addition, constant envelope modulation such as PSK has 
received particular interest recently with the emergence of 
large scale MIMO systems J35| . All the above motivate our 
focus on PSK constellations. Finally, we note that, in line with 
the relevant literature we assume a time division duplex (TDD) 
transmission here, where the base station directly estimates 
the downlink channel using uplink pilot symbols and uplink- 
downlink channel reciprocity. 

The rest of the paper is organized as follows. Section II 
introduces the system model considered in this paper and 
briefly describes the two conventional optimizations of inter¬ 
est: power minimization and SINR balancing. In section III 
the proposed beamforming optimization based on constructive 
interference is detailed for the case of power minimization, 
while its multicasting equivalent optimization is shown in 
section IV. Section V presents the constructive interference 


optimization for the SINR balancing problem. The CSI robust 
versions of these optimizations are examined in section VI for 
both power minimization and SINR balancing, for the case of 
bounded CSI errors. Finally numerical results are illustrated 
and discussed in section VII and concluding remarks are 
given in section VIII. In the following, we use the terms 
beamforming and precoding interchangeably, in-line with the 
most relevant literature. 


II. System Model and Beamforming Optimization 


Consider a K— user Gaussian broadcast channel where an 
7V-antenna BS transmits signals to K single-antenna users. 
Channel vector, precoding vector, received noise, data and 
SINR constraints for the i —th user are denoted as hj,G, 
rii, di = de^ i and T^ respectively with d denoting the 
constant amplitude of the PSK modulated symbols and rii ~ 
CJ\f(0, Wq), Vi where CA/’(/i,cr 2 ) denotes the circularly sym¬ 
metric complex Gaussian distribution with mean p and vari¬ 
ance a 2 . The received signal at user i is 

K 

Vi = h 1 k d k + Hi 
k =1 
K 

= hf^t k e j ^-^di+ni. (1) 

k=1 


The receive SINR at the i-th receiver for this scenario is 
given as 


\nfti 


H = 


J2k=i,k^i 

The transmit signal is 


+ AV 


Vi. 


( 2 ) 


K K 

yyt k dk = 

k =1 k =1 


( 3 ) 


where any of the users’ symbols di = de can be taken as 
reference. Without loss of generality, let us use d\ = de^ 1 
as reference hereof. Accordingly, the instantaneous transmit 
power for a signal with envelope d = 1 is defined as 

2 

. (4) 


Pt = 


k=l 


A. Power Minimization 

The conventional power minimization precoder, treating all 
interference as harmful, aims to minimize the average transmit 
power, subject to interference constraints as shown below fl6| 


min 

{ti} 

£ 11*4 

s.t. 



(5) 


\hfti 


Efc=i,fc^i + N 0 


>r i,v*. 


While optimal from a stochastic viewpoint, the above pre¬ 
coder ignores the fact that, instantaneously, interference can 
contribute constructively to the received signal power (9), and 
therefore from an instantaneous point of view we later show 
that it is suboptimal. 
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Fig. 1. Constructive interference in a) BPSK, b) QPSK and c) 8PSK 


B. SINR Balancing 

As regards the second optimization that is of interest in 
this work, SINR balancing maximizes the minimum achievable 
SINR subject to a transmit power budget, in the form 


max 

t i 

s.t. 


\tifU 


£lWI 2 <-P 




i =1 


( 6 ) 


where P denotes the total transmit power budget. We note that 
the above optimization is non-convex and the solution involves 
more complex iterative solutions ED- 


III. Proposed Optimization based on Constructive 
Interference and Phase Constraints 

With the knowledge of the downlink channel and all user’s 
data readily available at the transmitter, and with the aim 
of exploiting the instantaneous interference, the interference 
for PSK modulation can be classified to constructive and 
destructive based on known criteria (9), (32) . For clarity, these 
are summarized schematically in Fig. [T] for the basic PSK 
constellations. Here the scalar 7 denotes the threshold distance 
to the decision variables of the constellation, that relates to 
the SNR constraint, as detailed in the following. In brief, 
constructive interference is defined as the interference that 
moves the received symbols away from the decision thresholds 
of the constellation. This represents the green areas in the 
constellations of Fig. [T] where these are taken with reference 
to a minimum distance from the decision thresholds as per 
the SNR constraints. Note that these need not center on the 
nominal constellation points (the black dots in the figure) for 
generic SNR constraints. We refer the reader to 
for further details on this topic. 


As per the above classification and discussion, the opti¬ 
mizations in can be modified to take the constructive 

interference into account in the power minimization. This can 
be done by imposing interference constraints, not in terms of 
suppressing the stochastic interference, but rather optimizing 
instantaneous interference to contribute to the received signal 
power, thus reducing the required transmit power accordingly. 
For the case when interference has been aligned by means of 
precoding vectors tfc to overlap constructively with the signal 
of interest, all interference in 0 contributes constructively 
to the useful signal. Therefore all interference terms form 
components of the useful signal energy, which is given by 
the squared magnitude of the full sum in 0 . Accordingly, it 
has been shown in 0 that the receive SNR is given as 

hf 1 t>kdk 


where all interference contributes in the useful received sig¬ 
nal power. Accordingly, and based on the classification criteria 
detailed in [8] and Fig. [T] for M-PSK based on constructive 
interference, where the received interference is aligned to the 
symbols of interest, the problem can be reformulated as 


min 

ft*} 


s.t. 


K 


y7 fc e i(0,c ~ i;f>i) 

k =1 

Z(hf^t fc 4) =Z(di),Vi 


k=l 


Re > 7lVVo,Vi, (8) 


where Re(x), lm(x), sign(x) and Zx denote the real and 
imaginary part, the sign and angle of x respectively. Here the 
angle of interference is strictly constrained to equal the angle 
of the symbol of interest. The problem can be equivalently 
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Fig. 2. Optimization regions for (a) conventional precoding and (b) precoding for interference exploitation and generic optimization constraints (QPSK 
example). 


formulated as 


K 


min 

k=l 


( K 

s.t. 

Im ( hf yf 


\ k =l 


( K 


Re (hf 


We note the use of the sum of phase shifted (by the phase 
of the symbol of interest <pi) interfering symbols plus the 
symbol of interest in the above expressions. This serves to 
isolate the received amplitude and phase shift in the symbol of 
interest due to interference. Here, the first constraint requires 
that the interference perfectly aligns with the phase of the 
symbol of interest, to ensure that it overlaps constructively to 
the useful symbol © (32). The second constraint requires 
that this constructive interference is enough to satisfy the 
receive SNR threshold. Note that the above two conditions 
contain K equations and K inequalities, while there are 2 N 
real variables. Therefore, for the case N > K there are 
always sufficient degrees of freedom to satisfy these two sets 
of constraints, while our results in the following show that 
the proposed can be feasible with high probability even for 
cases with N < K. Still, it can be seen that due to the 
strict angle constraint, the formulation ([9]) is more constrained 
than the constructive interference regions in Fig. [T] where the 
strict angle constraints do not exist. To obtain a more relaxed 
optimization for M— PSK, let us observe the constellation 
example shown in the diagram of Fig. [2] for QPSK. Here, we 
have used the definitions olr = Re (hf J2k=i 

and Oil = Im (hf which are the real and 

imaginary components of iji = hf Ylk=i in the 

figure, the received symbol ignoring noise, phase shifted by 


the phase of the desired symbol. We also define 7 = \/rvVo. 
By means of their definition, aR and aj essentially shift the 
observation of received symbol onto the axis from the origin 
of the constellation diagram to the constellation symbol of 
interest. Clearly, aR provides a measure of the amplification 
of the received constellation point along the axis of the 
nominal constellation point due to constructive interference 
and ai provides a measure of the angle shift from the original 
constellation point, i.e. the deviation from the axis of the 
nominal constellation point with phase fa. 

At this point, we should emphasize the key idea in the 
proposed optimization which is the central strength of the 
proposed scheme that relaxes the optimization and allows 
additional reduction in the transmit power. In conventional 
optimizations, the signal power is optimized subject to SINR 
constraints. This is equivalent to constraining the interference 
each user experiences, so that the received symbol is within a 
certain distance from the nominal constellation symbol. From 
the view point of multiple users this essentially constrains 
the transmit vectors such that the received symbol y is con¬ 
tained within a circle (or a hyper-sphere for multidimensional 
optimizations) around the nominal constellation point d , so 
that the interference caused by the other symbols is limited. 
This is denoted by the dashed circle around the constellation 
point in Fig. [2ja). In contrast to this, here by use of the 
concept of constructive interference we allow a relaxation 
of aR and ai for all transmit symbols, under the condition 
that the interference caused is constructive, as secured by 
the constraints of the optimization. This gives rise to the 
constructive interference sector denoted by the green shaded 
sector in Fig. [2|b) 0, (32) . It can be seen that aR and aj 
are allowed to grow infinitely, as long as their ratio is such 
that the received symbol is contained within the constructive 
area of the constellation, i.e. the distances from the decision 
thresholds, as set by the SNR constraints Vi, are not violated. It 
is clear that the region in Fig. [2jb) is only constrained along 
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the vicinity of the decision thresholds, it therefore extends 
infinitely to the directions away from the decision thresholds 
and hence provides a more relaxed optimisation with respect 
to the conventional region of Fig. [2ja ). 

As regards the constructive area in the constellation, with 
respect to it can be seen that the angle of interference 
need not strictly align with the angle of the useful signal, as 
long as it falls within the constructive area of the constellation. 
For a given modulation order M the maximum angle shift in 
the constructive interference area is given by # = ±ir/M. 
Accordingly, to relax the optimization, ai is allowed to be 
non-zero as long as the resulting symbol lies within the 
constructive area of the constellation. Using basic geometry 
we arrive at the optimization problem expressed as in eq. ©• 
By noting that the last constraint actually incorporates the 
one for the real part of the received symbols, the optimization 
can be further reduced to the compact form of eq. ©• 
Clearly, the above optimization is equivalent to 


min 


s.t. 


K 




k=1 


,jO k-<f>i) 


K 


hf^tkdk - Z(di) 


k= 1 
K 


< $, \/i 


Re > 0VVo,V*- (12) 


k=1 


where for $ we have 


tanf? = tan# 1 — 


VTWo 


Re(hf £f=itfceW*-*«)' 


(13) 


It can be seen that the above optimization in © is more re¬ 
laxed than the zero-angle-shift optimization ([5]), which results 
in a smaller minimum in the transmit power. Moreover, it con¬ 
tains a number of K inequalities which result in an increased 
feasibility region compared to the conventional optimization, 
as will be shown in the following. Problem ( |TT| ) is a standard 
second-order cone program (SOCP), thus can be optimally 
solved using numerical software, such as Semudi. However, 
in the following section, we derive a more computationally 
efficient algorithm to solve it. For the illustrative example of 
BPSK used in the following results, the optimization can be 
modified to 


min 

ft*} 


s.t. 


K 


J2 tkt 

k =1 


J(<f>k—<f> l) 


K 


(14) 


R e hf y^t fe d fe sign(di) > vU/Vo, V«, 


k=1 


In ( p~4| ) the constraint requires that the interference on each 
user’s symbol has the same sign as the symbol of interest (and 
is therefore constructive) and that this constructive interference 
has enough power to satisfy the SNR threshold. For the case 
of QPSK the same principle needs to be applied separately to 


the real and imaginary part of the received symbol (see Fig. 
[TJb)) and hence we have 

2 


mm 


K 


tfe( 

k=1 


l) 


K 


(15) 


™ Vi. 


s.t. Re |^hf y^tfcdfcj sign(Re(dj)) > 

Im ^ tkdi^J sign(lm(di)) > / r ^° 


Vi. 


IV. A Virtual Multicast Formulation and a New 
Efficient Algorithm 


A. A Virtual Multicast Formulation of © 

The optimization in © can be re-cast as a virtual multi¬ 
casting problem (36]] according to the following theorem. 

Theorem 1: The broadcast problem ( pTj ) with constructive 
interference is equivalent to the multicast problem below 

min || w || 2 (16) 

w 


S.t. 



< (^Re — v/r 2 Ao) tan#, Vi, 


where the modified channel is defined as = h^e- 7 ^ 1 ( t >i \ 
To be more specific, the optimal solutions to © and ©> 
{t*} and w* respectively, have the following relation 


t* 

n 

t k 


w* 

~K ’ 

e j(4>i-4>k) 


w * e iOi -<f>k) 

K 


k = 2,- 


(17) 
,K.( 18) 


Proof: First we re-write the constraint in fTT] ) as ©. 

Observe that with ©, the composite precoding term 
Ef=i t k e^~^ in © can be treated as a single vector w 
precoder, i.e., Ylk=i t k e^^ k ~^ = w. Therefore the multicast 
reformulation © follows immediately. 

Suppose the optimal solutions of user l’s precoding vector 
in fTT) is t*. Without loss of optimality, the other users’ 
precoding vectors can be simple rotated versions of t*, i.e., 
t* = i = 2, • • • , K. As a special case, we have 

^ = tf This completes the proof. ■ 

Note that different from the classical multicast beamforming 
design, which is non-convex and difficult to solve (36j, © is 
a convex problem with a quadratic objective function and 2 K 
linear constraints thus easily solvable. In a similar fashion, the 
problem (l2| ) will have a similar multicast formulation. 

Theorem 1 provides useful insight into the structure of 
the precoding vectors. It tells us that the original broadcast 
problem now becomes a multicast problem if interference is 
utilized constructively. This is understandable because we take 
into account the correlation of the originally independent data 
streams, therefore the transmission can be re-formulated as 
sending a common data stream to all users, and re-shaping 
the channel and the resulting signal overlap, so that the 
intended data is delivered to each receiver. More importantly, 
the multicast problem contains only a single vector w and 
is therefore more easily solved than the original broadcast 
problem. The rest of this section is devoted to deriving an 
efficient algorithm to solve ©. 
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mm 

{t,} 


s.t. 


K 


53 tfee ' 




k = 1 


Re (f? 51 > \/iVVo,V* 


K 


Im h f 53 tfc< 


j(<f>k 4i) 


k=1 


iV 


< ^Re |^hf 


— \/rVVn I tan 0, Vi 


( 10 ) 


mm 

{ti} 


s.t. 


K 


53* fe e- 


, j (( j > k -< f > l ) 


k=l 


K 


Im h f 53 tfc< 


J(4>k~4>i) 


k=l 


K 


< Re hf - yfVYo tan 9, Vi 


k=l 


(ID 


K 


Im I hfe- 7 ^ 1 £t*e* 


(0fe-0l) 


fc=l 


K 


< Re V - y^Ao tan (9, Vi. 


/c=l 


(19) 


/?. Tfea/ Representation of the Problem 

For convenience, we separate the real and imaginary parts 
of each complex notation as follows 

hi = h Ri + jhn, w = w# + jw/, (20) 

where = Re(h i ),h/ i = Im^^Wi? = Re(w),wj = 
Im(w), j = V^T. 

Further define real-valued vectors 

fi = [h^-h/^wi = [w/;w fl ],w 2 = [w^; -w 7 ] (21) 

It is easy to verify that wi = IIw 2 , where II = 
[Oat In] -In Oat] - Oat and In denote N xN all-zero matrix 
and identity matrix, respectively. With the new notations, we 


mm 

{w 2 } 


w 2 


Wo = 


Ef =1 (^-^)gf-Ef =1 (^ + ^)ff 


It is not surprising to see that the optimal w 2 is the linear 
combination of the channel coefficients. 


Substituting (27) into £(w 2 ,/x, i/) leads to 


C(n,v) = 


\E?=M-Vi)gi-E“=M + Vi)V 


T ||2 


K 


+ ^/TiNo tan 6 53 (w 


(28) 


i = 1 


Define A = [/x; v\ and rearrange the terms, and we can 
obtain the following dual problem: 


can express the real and imaginary parts in (17) as follows 
Re(hfw) = f/ 1 w 2 ,Im(hfw) = f z T wi = f 2 T IIw 2 . (22) 


IIAAII 


max 

A>0 


v / TVVoi r A, 


(29) 


As a consequence, the constraint in (17) can be rewritten as 

|f 2 T IIw 2 1 < f J w 2 tan 6 — y/TiNo tan #, Vi. (23) 

Define = II T fi, and we rewrite (16) as 
|| 2 


s.t. gfw 2 — f 2 T w 2 tan# + y^Ao tan 0 < 0, Vi,(24) 
— gf w 2 — w 2 tan 0 + y^TVo tan 0 < 0, V(i25) 

C. The Dual Problem 

To further simplify the optimization let us formulate its 
dual problem. We let /i, v > 0 be the dual vector variables 
associated with the two sets of constraints in ( |24| ) and ( [25] ) 
respectively, and consider the Lagrangian of ( [24] ) as ( |26| ) 

The dual objective is thus given by min W2 £(w 2 , /x, v). 
Setting = 0 , we obtain the structure of the optimal 

w 2 below: 


where A = [f — g f + g] is an 2 A x 2 K real matrix, and the 
i-th column of f and g are defined in (21) and after (23). 

The problem ( [29] ) is a non-negative least-squares problem. It 
has wide applications but is known to be a challenging problem 
(38) . Without the non-negative constraint, its solution is given 
by 

A* = 2v / fAo(A T A)- 1 l. (30) 

Based on the above results, we have the following corollary 
to provide useful insight. 

Corollary 1: If A* is a positive vector , the optimal solution 
to the original problem ( [IT] ) with constructive interference 
is achieved when Im Thf = 0,Vi. The 


optimal precoding solution can be obtained from (30) and {21). 
This is can be explained by the fact that A is the dual variable 
and needs to satisfy the complementary slackness conditions 


in 


Mi(gf w 2 - ffw 2 tan<9+ V / TiAotan<9) = 0, Vi, (31) 
Vi{— gfw 2 - ffw 2 tan# + v/TVVotan#) = 0, V{32) 


2 


(27) 
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C(w 2 ,n,v) 


K 

IIW 2 II 2 + Mi (gf W2 - ff w 2 t an e - \/lVVo tan 9) (26) 

i= 1 
K 

+ za gf w 2 — f 2 T w 2 tan 9 — \/I\Nq tan 6^ 
i=1 


K 


K 


K 


w 2 || 2 + ( - ^i)gf - + ^) f f tan<9 J w 2 + 0WWE(« + ^) 

\2=1 i=l / i=l 

w 2 1| 2 + y^/ij ((gf - ff tan0)w 2 + \/rVVo tan0^ + ^ za (-(gf + ff tan(9)w 2 + v / IVVo tan^ 


2=1 


2=1 


When fii > 0, za >0, this implies that 

gfw 2 - ffw 2 tan (9 + V r i^o tan (9 = 

-gf w 2 - f, T w 2 tan 9 + a/TV^o tan 0 = 0, (33) 

0 


and it follows that 

Im ( h f Ek= i t fc e»'W*-*)) = 0, Vi 


g» r w 2 


or 


V. Constructive Interference Optimization for 
SINR Balancing 

The above concept of constructive interference can be 
applied to the SINR balancing problem of <0i. The problem 
for the case of constructive interference can be reformulated 
as 


D. A Gradient Projection Algorithm to Solve 

The general solution to ( [29] ) is difficult to derive. Here 
we propose a gradient projection algorithm to solve it. The 
gradient projection algorithm is the natural extension of the 
unconstrained steepest descent algorithm to bound constrained 


problems |37J. To this end, we first rewrite ( [29] ) as a standard 
convex problem: 


mm /(A) = MAL _ y/r T N 0 \, 


(34) 


max 

ft i,r t } 


s.t. 


K 


hj"^2t k d k - Z( di ) 


k =1 
K 


<i?,Vi 


Re > vT^Vo,Vz 


k =1 


A 


5 Z tfc< 

fc=l 


<f>l) 


< P 


(37) 


where T = [Ti, • • • , Tx] T . It is easy to verify that the 
gradient of /(A) is given by 

V/(A) = - vTUVol. (35) 

We then propose the following Algorithm 1 to solve ( [29] ). Once 
the optimal dual solution A op t is found, the original precoding 
solution can be obtained from ( [27] ). 


which is not a convex problem because is concave. 

However, this can be simply resolved by replacing it with a 
new variable, i.e., V t 2 = \fTt- Then we obtain an equivalent 
problem of ( [38] ) 

It can be seen that, as opposed to the conventional SINR 
balancing problem ([6]), which is non-convex, the formulation 
in ( [38] ) is convex and can be solved by standard convex 
optimization techniques. 


Algorithm 1 Efficient Gradient Projection Algorithm to solve 

m _ 

l: Input: h, d,r,Ao 
2: begin 

3: Initialize arbitrarily A^ 0 ^ > 0 . 

4: In the n— th iteration, update A: 


- y/T T N 0 l\, o) 
7 (36) 

where the step size a n can be chosen according to the 
Armijo rule or some other line search scheme. 

5: Go back to line 4 until convergence. 

6: end 

7: Output: A op t. 


kO) _ . 


(n-!) 


A t AA ( " 


A. Multicast Formulation of © 

In a similar fashion to the power minimization problem, the 
optimization in ( [37] ) can be re-cast as a multicasting problem 
according to the following theorem. 

Theorem 2: Problem ( [38] ) is equivalent to the multicast 
problem below: 


max V t 2 
{w,r t2 } 


Im(hfw) <(Re(hfw) 

2 < P 


w 


■ r t2 V -/Vo ) tan 0, Vi, 
(39) 


where h, = l . To be more specific, if the optimal 

solutions to (37) and (16) are {t*} and w*, respectively, then 
































max 

{ti,r t2 } 

s.t. 


r t2 


K 


im u E tfc< 


)j( 0 /e ) 


k=1 


< ^Re^hf^t fe e^-^j -r t2V /]Vojtan0,Vi 


K 


E** 

*!=i 


,j(<t>k-<t> 1 ) 


< P 
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(38) 


they have the following relation: 


t 


* 

l 


t/c 


W* 

if ’ 


w * e i( 0 i- 0 /e) 


,fc = 2, 


(40) 
, if. (41) 


Proof: The proof follows the one for Theorem [I] ■ 

The above is a standard SOCP problem which can be effi¬ 
ciently solved using known approaches. 


VI. Robust Power Minimization with bounded CSI 
errors 

A. Channel Error Model and Problem Formulation 

In this section, we study the robust procoding design when 
the CSI is imperfectly known. We model user i’s actual 
channel as 

hi = hi + e;,VP, (42) 

where hi denotes the CSI estimates known to the BS, and e* 
represents the CSI uncertainty within the spherical set lii = 

{e*:|NI 2 <€K 

We assume that the BS has no knowledge about except 
for its error bound 5? thus we take a worst-case approach 
for the transmit precoding design to guarantee the resulting 
solution is robust to all possible channel uncertainties within 
Hi. The specific robust design problem is to minimize the 
overall transmit power Pt for ensuring the users’ individual 
SINR constraints by optimizing the precoding vector {C}, i.e., 

min P T s.t. SINRi > T^ , Ve^ G Hi,Mi. ( 43 ) 

w v y 


B. Conventional Robust Precoding 

In the conventional multiuser MISO systems, the total trans¬ 
mit power is Pt = Ylf=t 111 < || 2 and robust design problem can 
be expressed as 

ElNI 2 

{ i =l 

INTI 2 

S.t. K 1 • *' -> r^Vej eUu'ik. (44) 

£|hft fe | 2 + 7V 0 

k = 1 
k^i 


Theorem 3: The robust beamforming problem ( [44] ) can be 
relaxed to the following semi-definite programming (SDP) 
problem 


K 


min > trace (T») 

{T^O, Si > 0 }^ 
k= 1 


S.t. 


where 


h*Qihf - 7i^o - Si 6f hJQi 
Qihf Q, + Sfl 


h o vfc, 

(45) 


M 


Qi — Tj — Tj Tn Vfc. 


k = 1 
k^i 


The problem ( [45] ) is convex and hence can be optimally solved. 
The resulting objective value of ( [45] ) provides a lower bound 
for the conventional power minimization. 

Remark: When the SDP relaxation is tight, or ( |?5| ) returns 
all rank-1 solution {TO, then the optimal solution {to} to 
solve ( |4?|) can be obtained by matrix decomposition such that 
T i = t^t J, Vi; otherwise, the required power in (44) is always 
higher than that in m 


C. Robust Precoding based on Constructive Interference 

Based on the multicast formulation of the power minimiza¬ 
tion problem([T6]), the worst-case robust design for the case of 
constructive interference becomes 


s.t. 


Im — ^Re — y / fyVo^ tan# < 0, 

V || ei || 2 <6f,Vi. (46) 


The constraint in ( [46] ) is intractable due to the infinite 
number of error vectors. Below we show how to tackle it using 
convex optimization techniques. For ease of composition, we 
omit the user index. 


First notice that robust constraint (46) can be guaranteed by 
the modified constraint below: 

max ^ ^|lm ^h T w^ | — ^Re ^h T w^ — y / fiVoj tan (9^ < 0 . 

(47) 

We again separate the real and imaginary parts of each 
complex notation as follows 


The robust precoding design is characterized by the follow¬ 
ing theorem pO) . 


h = + j h 7 

= h R + jh T +e R + je T . 


( 48 ) 
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Further define real-valued error vector and channel estima¬ 
tion vector 

e = [e#; e T ], f = [h R - h/]. (49) 


Apparently ||e|| 2 < S 2 . With the new notations, we can re¬ 
express the real and imaginary parts in ([47]) as follows 


Im(h T w) 

= h|jW/+ h/w fl + e|jW/+ efw_R, 



= f T wi + e r wj; 

(50) 

Re(h T w) 

= - h/w/ + e^w_R - efw/, 



= f T W2 + e T W2- 

(51) 


As a consequence, ( [47] ) can be rewritten as 

max |h T wi + e T wi I — (h T W 2 + e T W 2 ) tan6> 

||e|| 2 <<5 2 ' ' V ' 

-Fy/rTVo tan 6 < 0. (52) 

The above constraint is equivalent to the following two con¬ 
straints: 

max h T wi + e T wi — (h T W 2 + e T W 2 ) tan # 
l|e|| 2 <<5 2 

+^/TNo tan 6 < 0, (53) 

max liwi — e T wi — (h T W 2 + e T W 2 ) tan# 

||e|| 2 <<5 2 

+ y/rWotan# < 0, (54) 

whose robust formulations are given by 

h T wi — h T w 2 tan# + S ||wi — w 2 tan#|| 

Ay/rWotan# < 0, (55) 

—h T wi — h T W 2 tan# + S ||wi + W 2 tan#|| 

+ y / TWotan# < 0. (56) 


Finally we reach the following robust problem formulation 

min || w i|| 2 

Wi,W2 

s.t. Constraints ( [55] ) and ([56]), Vi, 

wi = nw 2 . (57) 

Note that ( [57] ) is standard SOCP problem therefore can be 
efficiently solved. After we obtain the optimal w^,W 2 , the 
robust solution w* can be determined using the relation in 


D. Robust SINR Balancing 

For completeness, we also study the robust SINR balancing 
given total BS power constraint. Similar procedures can be 
taken as the above to derive it, therefore we give the problem 
formulation directly below 

min r t2 

Wi ,W 2 ,r t 2 

s.t. Constraints ( [55] ) and ([56]), Vi, 

wi = riw 2 . 

||wi|| 2 <P t . (58) 

The above is a typical SOCP problem that can be solved using 
standard approaches. Suppose the optimal V t 2 is rj 2 , then the 
maximum SINR value becomes T^ 2 . 



Fig. 3. Transmit power vs. ft for conventional and Cl, K= 4, N=5 



r t (dB) 


Fig. 4. Transmit power vs. Ft for conventional and Cl, K= 4, N=4 


VII. Numerical Results 

This section presents numerical results based on Monte 
Carlo simulations of the proposed optimization techniques, 
termed as Cl in the following, and conventional precoding 
based on optimization for the frequency flat Rayleigh fading 
statistically uncorrelated multiuser MISO channel for both 
cases of perfect and erroneous CSI. Systems with BPSK, 
QPSK and 8 PSK modulation are considered while it is clear 
that the benefits of the proposed technique extend to higher 
order modulation. To focus on the proposed concept, we com¬ 
pare the proposed Cl optimization solely to the conventional 
optimization of fl 6 | , (3l| , while it is clear the the benefits of 
the proposed concept extend to modified optimization designs 
in the literature, by direct application of the constructive 
interference concept and the adaptation of the relevant QoS 
constraints. We use ‘TV x iC to denote a mutiuser MISO sys¬ 
tem with N transmit antennas and K single-antenna terminals. 
Unless otherwise specified, QPSK is the default modulation 
scheme. 

First, in Fig. [3] we compare the average transmit power with 
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Fig. 5. Transmit power vs. Ft for conventional and Cl, K= 4, N=3 


r=10dB 



Fig. 6. Feasibility probability vs. N for conventional and Cl, K- 4, Tt=10dB 


the conventional optimization and the proposed optimization 
problems ( p~4| ) and © in the 5 x 4 scenario for BPSK and 
QPSK, respectively. Power savings of up to 50% can be seen. 
It can also be observed that the relaxed optimization CD 
results in significant power gains compared to the strict angle 
constraints in & The same trend can be observed for the 4x4 
systems case of Fig. [4] where the conventional optimization 
results in a solution with increased transmit power, because 
less transmit antennas are available at the BS. The gains in 
this case are amplified for the proposed relaxed optimization. 
The transmit power is also shown in Fig. [5] for the 3x4 
scenario where the conventional precoder is inapplicable, as 
the optimization in 0 is infeasible for N < K. Remarkably, 
the proposed optimization problem is feasible in the 92.6% of 
the cases. The relaxed nature of the problem indeed leads to 
larger feasibility regions. To illustrate the extended feasibility 
region for the proposed optimization problems, Fig. [6] shows 
the feasibility probability of a K = 4 user system with respect 
to the number of transmit antennas. It can be seen that while 
the conventional optimization is only feasible for N > K, the 


5xK 



Fig. 7. Average execution time vs. K for Cl optimization with the BC 
formulation of 411) , MC formulation of and gradient projection solver 
in Algorithm 1, N=5, QPSK 


proposed can be feasible for lower values of N. 

The complexity of the proposed power minimization prob¬ 
lem is addressed in Fig. [7] for the system with N = 5. 
While analytical complexity expressions are hard to derive 
for optimization-based precoding, here the complexity of the 
solvers with the broadcast (BC) formulation of © multicast 
(MC) formulation of ( fib] ) and gradient projection solver in 
Algorithm 1 is shown in terms of the average execution time 
of the optimization algorithm, with increasing numbers of 
users. While the BC and MC show a similar computational 
complexity, the gradient projection solver in Algorithm 1 
offers a significant reduction in the involved complexity down 
to less than 15%, which further motivates the multicast for¬ 
mulation of the optimization problem. Still, we note that the 
proposed optimization needs to be performed on a symbol-by- 
symbol basis. Accordingly, the proposed may involve excess 
complexity compared to conventional precoding optimization 
especially for slow fading scenarios, where the convectional 
channel-dependent precoding may not require the optimization 
to be performed frequently. Accordingly, for the slow fading 
scenario where the channel is constant over a transmission 
frame and for the example of an LTE frame with 14 symbol 
time slots, the proposed precoding with the gradient projection 
solver translates to a doubling of complexity per frame (14 x 
15% = 210%) for the proposed scheme w.r.t. conventional 
precoding optimization. However, for the obtained transmit- 
power benefits as shown in our results, the end complexity 
increase is definitely worthwhile the performance benefits 
offered. In fact, in terms of the ultimate metric of power 
efficiency of the transmitter, we note that for an LTE base 
station the transmit power is measured in the order of around 
20Watts, while the power consumption of the DSP processing 
is typically of the order of hundreds of milliWatts. Since 
with the proposed precoding our results show a halving of 
the transmit power at roughly double the DSP power, the 
gains in the power efficiency for the proposed scheme w.r.t. 
conventional precoding are therefore significant. 
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5x4 r=20 dB 




Fig. 8. Achievable Ft vs. Transmit power budget for conventional and Cl, Fig- 10. Transmit power vs. error bound for conventional and Cl, K- 4, 
K= 4, N=5 N=4, r = 20 dB. 


S 2 = 0.0002 




Fig. 11. Transmit power vs. SINR for conventional and Cl, K= 4, N=4, 
Fig. 9. Achievable F t vs. Transmit power budget for conventional and Cl, £2 _ ^q-4 
K=4, N=4 


Figs. [8] and [9] compare the achievable SNR for the SINR 
balancing problems <0i and ( [37] ). In Fig. [8] the achievable T t 
is shown for the 5x4 MISO where an SNR gain of about 
2dB can be observed. The same trend is observed in Fig. [9] for 
the 5x4 MISO with an SNR gain of about 3dB. It is worth 
noting that for the case of constructive interference these SNR 
gains are due to the effect of constructive interference. 


Finally Figs. [TO] and ED compare the performance of the 
proposed CSI-robust Cl precoder with the conventional CSI- 


robust precoder of |30| for the 4x4 MISO. Fig. [TO] shows the 
obtained transmit power with increasing CSI error bounds S 
where it can be seen that for values in the region of 5 2 = 10 -3 
the transmit power increases significantly. On the contrary, the 
proposed optimization shows a modest increase in transmit 
power for high values of S thanks to the relaxed optimization 
obtained by exploiting constructive interference. This is also 
shown in Fig. [TT] where the average transmit power is shown 
with increasing SNR thresholds, for S 2 = 10 -4 . Again the 
proposed shows a constant loss of less than ldB compared to 


the perfect CSI case, while conventional precoding shows an 
increased sensitivity to the CSI errors. 

VIII. Conclusion 

A number of improved optimization-based precoding de¬ 
signs have been proposed for the multi-user MISO downlink 
channel. By taking advantage of the constructive interference, 
the proposed schemes deliver reduced transmit power for given 
QoS constraints, or equivalently increased minimum SNRs 
for a given transmit power budget. Both optimizations were 
adapted to virtual multicast formulations which were shown to 
be more efficiently solvable. The proposed concept was further 
extended to robust designs for imperfect CSI with bounded 
CSI errors. In all cases the proposed schemes provide superior 
performance over the conventional precoding schemes, which 
confirms the benefit of constructive interference when the 
precoders are fully optimized. 

The concept of constructive interference opens up new 
opportunities for future work in the optimization of precod¬ 
ing designs. Firstly, the extension to QAM constellation is 






















































































12 


non-trivial, and it needs a careful redesign of constructive 
interference sectors for determining the relevant optimization 
constraints. Secondly, the robust design in Section VI takes 
a conservative approach while the optimal robust precoder 
design is an important but challenging problem to be studied. 
Thirdly, as this work considers a single-cell system, and 
as multi-cell cooperation is made practical in systems such 
as cloud-RAN (39) , it is worth investigating the potential 
of the proposed scheme in multi-cell environments. Another 
promising application scenario is multi-beam satellite systems 
where both CSI and data for different beams are available 
at the gateway station where forward link beamforming is 
designed (40| . 
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